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The spectrum and amplitude of the stochastic background of relic gravitons produced in a bounc- 
ing universe is calculated. The matter content of the model consists of dust and radiation fluids, and 
the bounce occurs due to quantum cosmological effects when the universe approaches the classical 
singularity in the contracting phase. The resulting amplitude is very small and it cannot be observed 
by any present and near future gravitational wave detector. Hence, as in the ekpyrotic model, any 
observation of these relic gravitons will rule out this type of quantum cosmological bouncing model. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

Bouncing cosmological models [l[ are being widely in- 
vestigated because, besides solving the singularity prob- 
lem in cosmology by construction, they can also solve 
the horizon and flatness puzzles, and lead to an almost 
scale invariant spectrum of scalar perturbations if the 
contracting phase is dominated by dust at large scales 
1- 

One of the calculations that have been done was the 
evaluation of the spectral index of long wavelength ten- 
sor perturbations, Ut, and in most of bouncing models 
they were found to be also scale invariant. With respect 
to their amplitudes, specific models must be worked out. 
For instance, in the cyclic ekpyrotic scenario, they were 
evaluated and it was shown that the amplitudes are too 
small to be detected by present gravitational wave detec- 
tors 0. 

In this paper we will calculate the spectrum and am- 
plitude of relic gravitons in a different bouncing cos- 
mological model. It consists of a Fricdmann-Lcmaitrc- 
Robertson- Walker (FLRW) universe filled by dust and 
radiation, which is contracting classically. As it ap- 
proaches the singularity, quantum cosmological effects 
on the background, here described through the Wheeler- 
DeWitt equation interpreted along the lines of the Bohm- 
de Broglie quantum theory, avoids the singularity and 
ejects the universe to the expanding phase we are now 
experiencing. Note that our bouncing model is very con- 
servative: there is nothing else than dust and radiation, 
we are working with a 3 + 1 dimensional space-time, and 
we are performing a canonical quantization of the second 
order perturbed (with respect to the FLRW background) 
Einstcin-Hilbert action of general relativity, interpreted 
along the lines of a quantum theory appropriate to quan- 
tum cosmology, namely, one which does not need any 
external agent to the quantum system to give a meaning 
to the quantum calculations. 



Evolution of quantum perturbations (scalar, vector 
and tensor) on these quantum backgrounds can be de- 
scribed by simple equations, as it was demonstrated in 
Refs. 0, [3. With these equations, we were able to cal- 
culate the spectrum (analytically and numerically) and 
amplitude (numerically) of its stochastic background of 
relic gravitons. Although, the spectrum and amplitude 
differ considerably from the cyclic ekpyrotic scenario, the 
main conclusion remains the same, namely, that the am- 
plitudes are too small to be detected by any present and 
near future gravitational wave detectors. 

The paper is divided as follows: in the next section 
we review the main aspects of the quantum cosmologi- 
cal bouncing model on which the relic gravitons evolves, 
and we obtain the dynamical equations that the tensor 
perturbations which describe these relic gravitons must 
obey. In section III we derive the expression for the criti- 
cal fraction of the relic gravitons energy density from the 
tensor perturbations described in section II. In section 
IV we calculate the spectrum and amplitude of this en- 
ergy density and the graviton strain either analytically 
(for the spectrum) and numerically. We end up with the 
conclusions. 



II. THE STOCHASTIC BACKGROUND OF 
RELIC GRAVITONS IN A BOUNCING 
QUANTUM COSMOLOGICAL MODEL 



We start reviewing the key ideas concerning tensor per- 
turbations in a perfect fluid quantum cosmological model 
in the Bohm-de Broglie interpretation as pushed forward 
in Refs. 0,0. 
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A. Hamiltonian formalism for tensor perturbations 

We start with an Einstein-Hilbert action coupled to a 
perfect fluid described by the Schutz formalism Q : 



S = S GR + SWiid = ^-^pT J V^gRd^x + J ^/~ r gPA i x, 

(1) 

where £ p] = (8itG n /3) 1 / 2 is the Planck length in nat- 
ural units (H = c = 1), P is the perfect fluid pressure 
whose density p is given by the equation of state P = ujp, 
with ui = const. The metric g in Eq. (fTJ) is decomposed 
into a background Friedman-Lcmaitrc- Robertson- Walker 
(FLKW) metric and into a first-order tensor perturbation 



Wij, and is given in the Arnowitt-Dcser-Misner (ADM) 
formalism by 

ds 2 = N 2 (r) dr 2 - a 2 hys (r) ( 7y + Wij ) darW. (2) 

The background metric 7^ is related to the spacelike hy- 
persurfaces with constant curvature /C (/C = 0,±1 for 
flat, open and closed space respectively), and lowers and 
raises the indices of the tensor perturbation ujy, which 
is transverse and traceless (i.e. , w lJ V = and w\ = 0, 
where the bar indicates a covariant derivative with re- 
spect to 7). N(t) is the lapse function and defines the 
gauge, once fixed. 

The second-order Hamiltonian for the gravitational 
model described by the action ([T]) and metric ([2]) can 
be written as (see Ref. [f| for further details) 
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where the quantities P a , n iJ , P T are the momenta canon- 
ically conjugate to the scale factor, the tensor perturba- 
tions, and to the fluid degree of freedom, respectively. 
These quantities have been redefined in order to be di- 
mensionless. For instance, the physical scale factor a p h ys 
can be obtained from the dimensionless a present in (J3j) 
through a p h ys = t^afy/V, where V is the comoving vol- 
ume of the background spacelike hypersurfaces. This 
Hamiltonian, which is zero due to the constraint Hq w 0, 
yields the correct Einstein equations both at zeroth and 
first order in the perturbations, as can be checked explic- 
itly. In order to obtain its expression, no assumption has 
been made about the background dynamics, just Legen- 
dre and canonical transformations have been performed. 

In the quantum regime, this Hamiltonian can be sub- 
stantially simplified through the implementation of the 
quantum canonical transformation generated by 



Performing an inverse Legendre transform on the 
second-order piece of the Hamiltonian (jSJ, and restoring 
the constant £ rl , we get the following Lagrangian density 



£ (2) 



24£ 2 



2K.w l jW J 



(.0) 



where ^g a p is the background piece of the full metric 
([2]). Note that Lagrangian (|6|) coincides with the one 
derived in [|[ for classical backgrounds. 



U = cxp(zG q ) = exp 



12 1 



(4) 



where 



= \ (P a a + aP a ) and Q = J d : 
are the self-adjoint operators associated with the corre- 
sponding classical variables, yielding, for a particular fac- 
tor ordering of (see Section 3 in Ref. |6| for further 
details), the following quantum Hamiltonian 
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B. Quantum evolution of the background and 
perturbation variables 



The quantization procedure of the background and 
tensor perturbations can be implemented by imposing 
H ty(a,Wij) = 0. The Wheeler-DeWitt equation in this 
case reads 
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where we have chosen T as the time variable, which is 

Next, us- 



equivalent to impose the time gauge N = a 3 ", 
ing the Bohm-dc Broglic interpretation of quantum me- 
chanics @, making the separation ansatz for the wave 
functional ^[a,Wij,T] = y(q, T)ijj\a, wg, T], and follow- 
ing the reasoning of Ref. [10| ; we can show that Eq. (0 
can be split into two, namely 



Using the Bohm-de Broglie interpretation, Eq. (j8]) can 
now be solved as in Refs. [IlT[l3| . yielding a Bohmian 
quantum trajectory a(T), which in turn can be used to 
simplify Eq. ([9]) Q. Indeed, as one can understand a(T) 
as a prescribed function of time, one can perform the 
time dependent unitary transformation 



.dip a'^- 1 d 2 ip 
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and 
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U = cxp< i 
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yielding the following simple form for the Schrodinger equation for the perturbations: 
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(11) 



A transformation to conformal time rj, a 3uj ~ 1 dT — di], was also done, and a prime ' denotes the derivative with respect 
to r) (see again Ref. [f| for details). 

This Schrodinger equation is identical to the one used in semi-classical gravity for linear tensor perturbations, 
but here it was obtained without ever using the background Einstein's equations. Hence it can be used when the 
background is also quantized. It is important to stress that the function a{rf) which is present in Eq. (jlip is no 
longer the classical solution for the scale factor, but the quantum Bohmian solution. This fact leads to some different 
consequences with respect to the usual semi-classical approach, as we will see. 

In the Heisenberg representation, the equations for the operator Wij read 



^ + 2-^-^+2^=0, 



(12) 



which corresponds to the usual equation for quantum tensor perturbations in classical backgrounds (8| ■ 

We can proceed with the usual analysis, but now taking the quantum Bohmian solution 0(77) coming from Eq. © 
as the new pump field. In order to obtain these background quantum solutions, we follow the procedure discussed in 
Ref. [|[ which, for fiat spatial sections (K. = 0), yields the following guidance relation 



da 
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~~2 ~da' 



in accordance with the classical relations da/dT = {a,H} 
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dS/da. For the wave function 
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which comes from the initial normalized gaussian at T = 



*'"""<*> -{ni) -»l-5il- (15) 



where TJ, is an arbitrary constant, the Bohmian quantum 
trajectory for the scale factor is given by 



III. STOCHASTIC BACKGROUND OF RELIC 
GRAVITONS 



a(T) = a b 



(16) 



where a& is the minimum value for the scale factor at the 
bounce T = 0. Note that this solution has no singular- 
ities and tends to the classical solution when T — > ±oo. 
The quantity TJ, gives the curvature scale at the bounce, 

T = n^ w 

^bounce — b^fo 

The tensor perturbations are quantum-mechanical op- 
erators, hence it is convenient to expand them into 
Fourier modes and subject them to quantization rules: 

d 3 k 
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*(A) 



E 

A=+,x 
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where x = (77, x), e-^' = e-^' ^kj is the polarization ten- 
sor for the two graviton polarization states + and x la- 
beled by A, and satisfies 



e Wij (V) = 2S 



XX' 



(18) 



Also, (77) are mode functions, and a k A ' ' , a k ' ' are cre- 
ation and annihilation operators, respectively. Such op- 
erators satisfy the equal-time commutation relations 



'k 



~(A) .(A')t' 
ffl k .°k' 



= 5 XX ,S^ ( k - k') 



(A) ~(A') 



k' 



(A)t ~(V)t 



and the quantum vacuum is defined by 



4 A) |o> = 0. 



Next, inserting the Fourier expansion into Eq. 
we get the mode equation 



(19) 
(20) 

(21) 
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Introducing the canonical amplitude as 
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,(A) 



mode equation (|22|) becomes 
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2K-^ 
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0. 



(22) 



(23) 



(24) 



for each graviton polarization state. In the present work, 
we will only be concerned with the case K, = 0, i.e., a 
Universe with flat spatial sections. 



In the ACDM model, quantum fluctuations arising 
during inflation lead to a nearly scale-invariant spectrum 
of density (scalar) [l4| and gravitational waves (tensor) 
perturbations [l5|. As in the case of CMB, we also expect 
relic gravitons generated in this early epoch to form a 
background - the stochastic background of relic gravitons 
(SBRG) - to be hopefully detected by the high-sensitivity 
gravitational waves (GW) detectors. The physical ob- 
servable to be measured by such GW detectors is the 
critical fraction of the relic gravitons energy density Pgw 
given by (see [lfl [I?} and references therein) 



GW 



(r?) 



Pgw (rf) 



(25) 



where p c = (H/£ m ) is the critical energy density. The 
energy density paw carried by the relic gravitons is sim- 
ply the 0-0 component of the stress energy-momentum 
tensor of the tensor perturbations, that is Pgw — T , 
where 



T. 



a/3 



(26) 



and £( 2 > is the Lagrangian ^ for the tensor perturba- 
tions (recall that Lagrangian (|6|) holds either in classical 
or quantum backgrounds). Considering only spatially- 
flat models (K, = 0), the energy density of the relic gravi- 
tons can be derived by inserting © into (|2"6"|) and taking 
the 0-0 component, which yields the classical expression 



Pgw 



12£ 2 



;W J Wi 



(27) 



Since tensor fluctuations generated during inflation are 
of quantum nature, the energy density of the relic gravi- 
tons becomes the expected value of the operator T° in 
the vacuum state defined in (|2"TT): 



Pgw 



(0|T° |0). 



(28) 



Substituting the Fourier expansion (fTTj) into (f2"7} . and the 
result into (|28|) . we find, using the commutation relations 
(119I20I21I) . that 
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In terms of the canonical field v ( k ' introduced in (|23l) , where we have introduced the energy density per mode 



expression 



turns into 



Pgw (rj) = dink pgw (k, 77) 
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where % = a' /a is the Hubble parameter in conformal 
time coordinates. 

From pil) the energy density parameter (|25p per mode 
is given by 



&gw (k, rf) 



2it 2 a 2 H 2 
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(32) 



By solving the mode equation (|24p and inserting the 
result into (|32[) . we promptly determine the energy den- 
sity parameter per logarithmic frequency interval v = 
k/(2ira), which is the physical quantity to be confronted 
with future observations. 



IV. STOCHASTIC BACKGROUND OF RELIC 
GRAVITONS IN A DUST-RADIATION 
BOUNCING UNIVERSE 

A. Scale Factor 

The scale factor for the radiation dominated bounce 
will be given by Eq. (fTo) , with w = 1/3. In this case, we 
have dT = dr] so, we find 



dbou 




(33) 



where r\b is a free parameter that determines the duration 
of the bouncing phase. 

The scale factor before and after the bounce is the 
usual classical solution for a Universe filled with radiation 
and dust [H, 



O-radWl = a e q 



T2- 



V* 



(34) 



where the minus and plus signs refers to the epochs before 
and after the bounce, respectively, a eq is the value of the 
scale factor at matter-radiation equality, and 



V* = 2Rh 



l + z, 



(35) 



where Rr = ^-/{clqHq) is the co-moving Hubble radius, 
do and Ho are the present values of the scale factor and 
Hubble radius, respectively, and z eq ~ 3 x 10 3 is the 
value of the redshift at matter-radiation equality. As the 
bounce should occur much before the matter-radiation 
equality, then rjb << 

The above expressions for ab ourace and a ra d can be con- 
densed into 



a(77) = a eq 



Vb 



(36) 



which promptly recovers Eqs. (f3"3"|) and (f3"4"]) in the limits 
77 << ?7* and 77 >> 77^, respectively, and assuming that 
cib = 2a eq rjb/i]*. 

With this expression for the scale factor, we are ready 
to solve the mode equation, Eq. 



B. Mode Equation 

For numerical purposes, we prefer to work with the 
normalized variables 



V 



77 
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2R H 
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In terms of these variables, the mode equation for a 
flat Universe becomes 



d 2 V n 

dfj 2 



da 1 
dfj 2 a 
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where we are treating separately each possible graviton 
polarization state, so we have dropped the superscript 
(A). 



6 



The scale factor given in Eq. (l36l) becomes 



a{r]) = a eq 



1 



where rji, = t]i > /2Rh is a free parameter and 
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0.018. 
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In the original set of variables rj and A:, the initial con- 
ditions to solve the mode equation are the vacuum con- 
ditions to be set at r] ini = — oo: 
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For fj = fjo, this expression reduces to 
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Using the new variables fj and n, the initial conditions 
to solve Eq. (|38|) for one of the two polarization modes 
are 



Rh _ 
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dVn 
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C. Analytical Results 

One can anticipate the spectral behavior of the above 
quantities through some analytical reasonings. Note first 
that the spectral index of tensor perturbations at the 
moment they get smaller than the curvature scale in the 
expanding phase of the above quantum bouncing models 
was calculated in Ref. Nj, and it reads 



riT 



12w 
l + 3w y 



(49) 



where w is the equation of state parameter of the fluid 
which dominates the cosmic evolution when the mode is 
getting bigger than the curvature scale of the Universe in 
the contracting phase. When the fluid is dust, tit = 0, 
when it is radiation, = 2. 
The power spectrum is 



We have checked that we can, in this case, numerically 
consider — oo as any value of < — 3. 

The present value of fj, fjo, is obtained by solving Eq. 
(f3"9"| for fj = fjo, which results in 



1 + z, 



0.982. 



(45) 



We can constrain the free parameter fjb by demanding 
that the end of the bouncing phase happens between the 
Planck time [z ps 10 32 and the beginning of nucleosyn- 
thesis (z ps 10 10 ), which results in 



10~ 31 < fjb < 10" 9 . 



(46) 



We have solved Eq. (|38|) from r)i n i to r/o for four dif- 
ferent values of fjb in this range: 10~ 30 , 10 -24 , 10~ 18 and 

io- 12 . 

In terms of the new variables fj and n, the expression 
for the density parameter in Eq. (|32p , for one of the two 
possible polarizations, can be rewritten as 



V oc n 3 |/i„| 2 = (Ah n ) 



2 oc n nT , 



(50) 



where Ah n is the strain of the gravitational wave with 
wavenumber n. Hence, one can obtain the spectra of 
the strain at the moment rj c of curvature scale cross- 
ing at dust domination and radiation domination from 
Eq. (|49[) , and they are Ah n (r] c ) oc n° and Ah n (r) c ) oc n, 
respectively. 

However, we need the strain today, Ah n (rjo). Note 
that the curvature scale is given by l c = (a 3 /a") 1 / 2 
(which is equal to the Hubble radius in the case of a 
perfect fluid dominated classical model) , and it is compa- 
rable to the physical wavelength of the mode k = n/ Rh, 
^phys = a/k, when k 2 = a" /a. Hence, for Z p h ys < lc, or 
k > a" I a, Vk oscillates (see the mode equation (|2~4"|) 
for K. = 0). Therefore, as Vk = a^k, one has that 
lMfc|( 7 7o)a( 7 7o) = lMfc|( 7 7c)a(77 c ) after curvature scale cross- 
ing at the expanding phase. Note that a(rj c ) as a function 
of k can be obtained through 



lphys{ih) = a(rj c )/k = l c (rj c ) oc Rh(Vc) oc t(r) c ) oc a 3/2 (r) c ) 



(51) 
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when rj c is in the dust dominated era, and 

IphysiVc) = a(r] c )/k = l c (r) c ) oc R H {Vc) oc t(r) c ) oc a 2 {jj c ) 

(52) 

when rj c is in the radiation dominated era. Hence one 
gets that a(r] c ) oc k~ 2 and a(rj c ) oc fc _1 for dust and radi- 
ation dominated eras, respectively. Therefore, in terms 
of the frequency / = k/(2nao) = nHo/(4n), the transfer 
function for the strain, 



fW(/,r7o) <xf{Ah f (fj )] 2 . (57) 

These results are confirmed by the numerical calcula- 
tion we show in the sequel. 



D. Numerical Results 



Ahf(r) )=T f Ahf( Vc ), (53) 
which is the same as the one for \fJ,f\, see Eq. (|50p . reads 



't 1 



fcq 



(54) 



where / cq = H cq a cq /a ~ 1CT 16 Hz and H ~ 1CT 18 Hz 
Finally, the potential V(fj) = (d 2 a/dfj 2 )/a of Eq. ([38l 
derived from Eq. ([3Uf reads 



F(?7) 



(55) 



and it has a maximum value at f/ = 0, Kn ax ~ % 2 - We 
see therefore that modes with n > I /fjb do not feel the 
potential at any time. Hence they are always oscillating 
with spectrum given by the initial vacuum state, 



v n (v) = 



(56) 



From Eq. (|50p . one then infers that the strain Ah n has 
spectrum n for n > 1/fjb, independently of r/. As the 
bounce itself has a duration much smaller then the clas- 
sical dust-radiation evolution, almost all modes that do 
not enter the curvature scale after the beginning of clas- 
sical radiation domination in the expanding phase have 
the above spectrum. 

Noting that our model is symmetric with respect to 
the bounce, implying that wavelengths which leave the 
curvature scale before the bounce at some fluid domi- 
nation era enter the curvature scale after the bounce at 
the same fluid domination, we can summarize the results 
concerning the spectrum of the strain in the following 
way: 

i) H ~ lCT 18 Hz < f < l(T 16 Hz K / cq ; 
Ahf(fjo) ex f- 2 f = r 2 . 

ii) f cq « 10- 16 Hz < / < / max = H /fj b ; 
Ah f (fj ) ck r 1 / = /°. 

hi) /max = H /r) b < /; Ahf(r) ) oc /. 

From the above spectrum, one can easily obtain the 
spectrum for the energy density parameter of gravita- 
tional waves from the relation 



Fig. [T] shows the value of the density parameter as a 
function of frequency / for the four different values of fjb 
we are considering. Fig. [2] shows the results of the strain 
amplitude 



[A/i„(?7o)] 2 = 6 ! 



(58) 



Both figures are in accordance with our previous ana- 
lytical considerations. 
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FIG. 1. Energy density parameter as a function of frequency 
for the primordial gravitational waves produced in our model. 
Each curve shows the results obtained for the indicated value 
of fjb. The dashed lines show for comparison two power laws 
proportional to n 2 and n 4 . 



In order to establish the possibility of detecting the 
primordial gravitational waves produced in our theoret- 
ical model, we show in Fig. |3] a zoom of the upper part 
of Fig. [I] together with other theoretical predictions and 
experimental sensitivities. 

We see that only for the lowest value of fjb we stud- 
ied, fjb = 10 -30 , and for a very high value of frequency, 
does our result reach the upper limits predicted by the 
inflationary model. Not even in this extreme case, is 
the stochastic background strong enough to be directly 
detected by present and future planned experiments. 
Nonetheless, we must emphasize that, while our calcu- 
lations are exact, the predictions presented for inflation 
are only upper limits and could, therefore, be much lower 
than the curve presented in Fig. [3] 
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FIG. 2. Strain amplitude as a function of frequency for the 
same primordial gravitational waves produced in our model. 
Each curve shows the results obtained for the indicated value 
of 776 . The dashed lines show for comparison two curves pro- 
portional to n~ 2 and n. 
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FIG. 3. Comparison of our results (blue curves at the bottom 
right) with experimental sensitivities (red curves) and a pre- 
diction the upper limits on the spectrum of primordial gravi- 
tational waves generated in inflationary models (black dashed 
curve). The red curves show the sensitivities achieved by 
LIGO's 5th run and the ones predicted for Advanced LIGO, 
the Einstein Telescope and LISA [llj. See and references 
therein. 



V. CONCLUSIONS 

In this paper we have calculated the amplitude and 
spectrum of energy density and strain of relic gravitons 
in a quantum bouncing cosmological model. The strain 
spectrum of this quantum bouncing model is different 
from the cyclic and inflationary scenarios. While these 
two models have spectra pa fc~ 2 ,« k~ x at dust domina- 
tion, and « fc ,« k° at radiation domination, respec- 
tively, our model have spectra « fc~ 2 and « k° at the 
same eras. 

One possible different scenario in which the power of 
relic gravitons could be enhanced could be obtained by 
adding to the matter content of the model some amount 
of stiff matter, which has primordial spectral index tit = 
2. This will be the subject of future investigations. 

As a final point, as in the cyclic ekpyrotic model, the 
resulting amplitude is too small to be detected by any 
gravitational wave detector. In particular, the sensitiv- 
ity of the future third generation of gravitational wave 
detectors, as for example the Einstein Telescope, could 
reach f^cw ~ 10~ 12 at the frequency range 10 — 100 Hz 
to an observation time of ~ 5 years and with a signal- 
to-noise ratio (S/N) ~ 3. Therefore, any detection of 
relic gravitons, in this frequency range, will rule out this 
type of quantum bouncing model as a viable cosmological 
model of the primordial universe. 
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